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SUMMARY 

A theoretical investigation has been made of the pressure distributions on and drag 
characteristics of the faces of forward-facing steps immersed in turbulent boundary layers at 
supersonic speeds. An approximate solution technique proposed by Uebelhack has been 
modified and extended to obtain a more consistent' numerical procedure. 

Numerical results obtained indicate that decreases in the assumed separation angle and 
increases in the assumed shear-layer spreading parameter result in decreased pressure levels 
and drag forces. Increases m the index of the assumed power-law velocity profiles caused 
small increases in the dividing streamline pressures at the smaller ratios of step height to 
boundary-layer thickness but had no significant effect on the drag. Modifications to the 
width of the shear layer and reverse-flow pressure distnbutions resulted in reasonably good 
agreement between theoretical and experimental face pressure distributions for moderate ratios 
of step height to boundary-layer thickness (from 2 to 7). By coupling a decrease in assumed 
separation angle with an increase in the shear-layer spreading parameter, the modified theory 
can be extended empirically to give good agreement between theory and experiment for 
ratios of step height to boundary-layer thickness as low as 0.05 at free-stream Mach numbers 
of 1.61 and 2.20. There is some disagreement between the results obtained m this investi- 
gation and those of Uebelhack which apparently cannot be ascribed to different methods of 
computation. 


INTRODUCTION 

The study of supersonic turbulent boundary-layer separation due to a forward-facing 
step has many useful applications. At low ratios of step height to boundary-layer thickness, 
the data and analytical procedures are helpful in assessing the effects of manufacturing 
discontinuities on aircraft drag and performance, at larger ratios, the results can be used to 
predict the effectiveness of spoilers and deceleration devices. The separation phenomena 
have been extensively investigated both theoretically and experimentally, yet in neither case 
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has much emphasis been focused on the flow at or near the step face From the standpoint 
of theory this lack of emphasis stems primarily from the complexity of the problem. The 
lack of experimental attention results from the very small models that usually must be used 
which makes it very difficult to install a sufficient number of static orifices to reasonably 
determine the pressure distribution over the step face. Inasmuch as the flow near the step 
face determines the drag characteristics of forward-facing steps, a need exists to investigate 
this flow region. A logical first step appeared to be a study of the utility of the theory 
proposed by Uebelhack (ref. 1), the results of this study are the subject of the present paper. 

Uebelhack’s theory is based on a highly simplified phenomenological model but appears 
to be in good agreement with experimental results. The theory is stated to be applicable 
only for ratios of step height to boundary-layer total thickness of greater than 2 (ref. 1) 

Also the theory provides an average pressure for the step face and a peak or reattachment 
pressure which is the stagnation pressure on the streamline dividing the flow passing over the 
step from that being recirculated in the separated-flow region In this study the implications 
of the theory as to the pressure distributions on the step face are examined and the require- 
ments for extending the usefulness of the theory to lower ratios of step height to boundary- 
layer thickness are established. In addition, variations m the basic parameters are investigated 
in more detail than by Uebelhack and some errors m procedure and experimental data inter- 
pretation are corrected. Calculations were made over a Mach number range from 1.61 
to 6.0, over a range of ratios of step height to boundary-layer total thickness from 0.001 
to oo, and over a range of boundary-layer velocity profiles from the 1/7 to the 1/11 power 
laws. 


SYMBOLS 

^0’^r^2’^3 arbitrary constants (eq. (13)) 

“a “2 e 

C„ Crocco number, = — — 

^ u u 

^max max 

Cj) drag coefficient 

d dividing streamline 

h step height 

j inviscid jet boundary streamline 

M Mach number 
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m 


mass, center of mixing 


Npj. Prandtl number 

n power-law exponent for velocity profile 

p pressure 

Rgj Reynolds number 

Tq total temperature 

u velocity along x-axis 

V velocity along y-axis 

X longitudinal coordinate in original coordinate system (fig. 20) 

Xq correction for mixing length (eq. (A25)) 

X* length of the mixing region (fig 1) 

X longitudinal coordinate in modified coordinate system or effective length of 

mixing region (fig 20) 

y lateral or step coordinate 

P shock-wave angle, degrees 

y specific-heat ratio 

6 boundary-layer thickness 

5* boundary-layer displacement thickness 

r? dimensionless coordinate, ya/x 

r) dimensionless coordinate, ya/x 
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d 


boundary-layer momentum thickness 


separation angle, degrees 
p density 

0 shear-spreading parameter (eq. (2)) 

^2 ^2 

0 velocity ratio, — = 

Subscnpts: 

a region of higher velocity in two mixing streams (identical to region above the 

shear layer behind the oblique separation shock) 

b region of lower velocity in two mixing streams, also region of reverse flow 

d dividing streamline 

e boundary-layer edge 

j mviscid jet boundary streamline values 

m center of mixing 

max maximum 

0 stagnation conditions 

ra,rb large positive and negative reference (fig. 18) 

1 ahead of oblique separation shock (free stream) 

2 behind oblique separation shock 
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THEORETICAL APPROACH 


The calculation of the pressures on the front of a forward-facing step is based on 
Uebelhack’s shear-flow model which is shown schematically in figure 1. As the flow 
approaches the step face, the boundary-layer flow separates from the plane wall, inducing an 
oblique shock. Behind the shock, the boundary layer is assumed to change to a free shear 
flow inclined at a constant separation angle 0^ to the surface. The shear flow expands 
as it proceeds downstream and reattaches to the step face with part of the flow passing 
over the step and part of the flow being reversed in direction in the separated-flow region 
near the wall. No physical mechanism is provided for reversing the flow. For analysis, the 
basic flow model is divided into the following four regions 

(1) Separation region 

(2) Free-mixing region 

(3) Reattachment region 

(4) Reverse-flow region 

The concepts and most important assumptions used to analyze the individual regions 
are described briefly in this section. (For a more detailed descnption of the method, see 
ref. 1.) A concise derivation of the basic equations is presented in appendix A and a 
description of the calculation procedure is presented in appendix B. Uebelhack’s derivation 
of the basic equations is repeated in this section since some modifications were made to the 
procedure as well as corrections of a number of typographical errors in reference 1. 

(1) Separation Region 

The flow separates from the wall before reaching the step face at a distance deter- 
mined by the separation angle 0^ (fig. 1). This angle is determined from correlation of 
experimental results (ref. 1). At the separation point an oblique shock is induced. Ahead 
of the shock, the boundary layer is characterized in terms of 5j, 5^, and 6y Behind 

the shock, the new boundary-layer characteristics are calculated by use of the oblique-shock 
theory. The boundary-layer momentum thickness behind the shock ^2 is the important 
parameter connecting mass and momentum flux. In terms of the initial boundary-layer 
characteristics and flow-separation parameter ^2 found to be 
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^2 
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cos 6^ 


sin 1 

fe - e] 

y SI 

1 cos 1 

- «s) 

' Sin j3 

cos j3 


- ( 6 , - 6 |)- 


(*^ - 
Sin /3 


cos 0g + 5j 


Pi - P2 

■^P2“le J 


where |3 is the oblique shock angle (fig. 1), Pj and P 2 are the static pressures ahead 

of and behind the shock, respectively, and M, is the Mach number outside the boundary 

z,e 

layer behind the shock 


(2) Free-Mixing Region 

The boundary layer behind the oblique shock is assumed to change instantaneously to 
a shear layer having the same momentum thickness as the boundary layer immediately behind 
the shock. For simplicity, the following assumptions are made. 

(1) The shear layer expands linearly. 

(2) The mixing occurs at the constant static pressure P 2 - 

(3) The shear-layer velocity distribution can be approximated by an error-function profile. 

(4) The flow is isoenergetic with unit Prandtl number ^Tq = Constant, Npj. = ij. 

The linearly spreading shear layer requires the use of a spreading or mixing param- 
eter a which is based on experimental correlations. For this investigation a was 
taken according to Korst and Chow’s suggestion for compressible fluids (ref. 2): 

a = 12 + 2.76M. (2) 

2,e 

Reference 3 presents a more recent correlation which yields values of a nearly twice as 
large as those indicated by equation (2). The recommended correlation is, however, highly 
tentative and intended for fully developed shear layers, whereas for practical ratios of step 
height to boundary-layer thickness (h/5 < 5), the experimental shear layer generally cannot 
be expected to be fully developed. Equation (2) yields results that appear to be in reason- 
able agreement with the still developing shear layers of reference 3. For h/Sj > 5, the 
shear layer may become more fully developed and higher a values may be appropriate. 
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The shear layer must be properly located’ specification that the momentum thickness 
of the shear layer be equal to that of the boundary layer places the shear layer at the 
correct longitudinal location; specification of the angle alines the layer at the desired 

angle with reference to the wall. Proper location of the center of the shear layer normal 
to the wall at the separation point requires two adjustments. The first of these is concerned 
with the fact that as the shear layer grows in width in the downstream direction from the 
virtual origin, the center of the layer, which is specified in terms of the symmetry of the 
velocity profile, is continuously displaced m a normal direction away from the initial position 
at the virtual origin. The adjustment is calculated by invoking the laws of conservation of 
mass and momentum within a control volume enclosing the upper and lower bounds of the 
shear layer. (See appendix A.) The second adjustment must account for the thickness of 
the boundary layer immediately behind the separation shock. This adjustment is accomplished 
by the assumed injection of air into the region below the shear layer to account for the 
momentum thickness of the boundary layer (See appendix A.) 

(3) Reattachment Region 

As the shear flow approaches the step, a part of the flow reattaches to the step face, 
while the remainder passes over the step (fig. 1) The streamline that identifies these two 
parts of the flow is defined as the dividing or reattachment streamline and the assumption 
IS made that this streamline reattaches at the outer comer of the step. The present theory 
does not provide any physical mechamsm for reversing the flow. In effect, the assumption 
is made that the flow in the reattaching shear layer below the dividing streamline can be 
stagnated without any turning of the flow or spreading of the streamlines. This procedure 
allows the retarding force on the step face to be calculated at the pnce of the unaccounted- 
for disappearance of the mass flow contained in this part of the shear layer. 

(4) Reverse-Flow Region 

In the Uebelhack separated-flow model, the flow beneath the shear layer is assumed to 
be reversed at constant velocity (fig. 1) with constant static pressure P 2 behind the sepa- 
ration shock. The mass flow m this region is equated to the mass flow in the reattaching 
shear layer below the dividing streamline at the step face. This procedure allows evaluation 
of the force on the step face adjacent to this region required to accelerate the subject mass 
flow from zero velocity to the final reverse-flow velocity; again the mass flow must appear 
unaccountably at the step face. 

Inasmuch as the experimental data do not indicate uniform reverse-flow velocities, the 
theory was modified for some calculations with the assumption that the pressure distnbutions 
m the reverse-flow region are mirror images of those in the shear-layer portion of the step 
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face with the focal or turning point located at half the step height Ihis approach will be 
defined as the “variable Tjj-b ” rnethod. 


final lk|uations 

One of the final results of the analysis is an integral over the pressure distribution on 
the step face. This integral is 


I 


1 p - p 


where 
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( 3 ) 


^d 

^rb 


longitudinal length of shear layer 
dimensionless coordinate for dividing streamline 

arbitrary dimensionless coordinate at which the error-function velocity is very 
small 


V 

h 


~ ^2 e/’^max’ isoenergetic Crocco number 

~ “ 2/^2 e’ velocity ratio 

= yo/x, dimensionless lateral coordinate 
height of the step face 


This pressure integral can be related to the free-stream static pressure existing ahead of the 
separation shock by 


rl ^ ^ ^ / ri P ^ ^ + A ^2 

Jo Pi ^ Vo P2 h y Pi 


( 4 ) 
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by 


From the preceding pressure integral, the drag coefficient on the step face is calculated 



(5) 


Uebelhack (ref. 1) did not consider the effective pressure distribution on the step face 
except for the stagnation pressure on the dividing streamline. To calculate this pressure he 
used the equation 


Po4 

Pi 



Pi 


( 6 ) 


where p^ ^ is the stagnation pressure at the end of the dividing streamline and Mj is 
the Mach number on the dividing streamline. The use of this equation, which applies to 
subsonic isentropic flows, is not compatible with the assumption of constant pressure and 
total temperature in the shear layer. Uebelhack does not specify his exact procedure for 
the case where the Mach number on the dividing streamline is supersonic. Because the 
implied pressure distributions on the step face were of interest in this investigation the 
correct equations for the pressures were derived from the pressure integral (eq. (3)) and are 
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for the shear-layer portion of the flow and 





for the reverse-flow portion of the step. Equation (4) is modified in equations (7) and (8) 
to the form 


Pi 



(9) 
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When the free-stream Mach number and step height are sufficiently large, a part of the 
shear layer in the separated-flow region is supersonic ^for example, for h^5j > 2 at 
Mj ^ = 5^ For all reattachment pressures m this case, but only for a few drag and pressure 
distnbution calculations, it was assumed that this flow was decelerated through a normal 
shock. The equations given in reference 3 were utilized to obtain the pressure ratio 
Pod/Pj when accounting for the normal shock. In the reverse-flow region, shocks cannot 
exist because the flow is accelerating. 

RESULTS AND DISCUSSION 


Constant 

Comparison with Uebelhack’s results.- A companson of the results computed in this 
investigation with those of Uebelhack is presented in figure 2. Here, the dividing-streamline 
stagnation or reattachment pressure p^ j and drag coefficient Cjj are plotted as a 
function of free-stream Mach number M, and the ratio of step height to free-stream total 
boundary-layer thickness (step-height ratio) h/6j. For a direct companson with Uebelhack’s 
results, the same basic approach with constant was utilized and the dividing-streamline 

pressure ratio was computed by the same isentropic flow procedure used by Uebelhack in 
his calculations. If the dividing-streamline Mach number is higher than 1, a normal shock 
is taken into account for the Pod/P i calculation. The calculations were made for a range 
of g from 2 to 6 and for a range of h/5j from 0.1 to 100 (essentially <»). 

The curves show that the present results are in only fair agreement with Uebelhack’s 
results. The Pod/P i curves of the present investigation are more linear and more 
uniformly spaced m both the Mach number and h/6j plots. In the case of Cp, 
Uebelhack’s results show more variation with h/6 j than the present calculations and, in 
fact, show an increase in Cp with a decrease in h/5j which is contrary to expectation. 
The reasons for the major part of these discrepancies, which are unexpected because the 
same basic equations are used for the calculations, are not known. 


Note that the curves indicate that Cp and Pod/P i asymptotically approach limiting 
values as h/6 ^ The reason for this is that for very large step-height ratios, the 

boundary-layer thickness is only a minute part of the step height. The boundary layer 
provides for the separation of the flow but the detail characteristics of the initial boundary 
layer fade into insignificance. Thus, Reynolds number effects may be expected to exist for 
low values of h/6j but little or no effects are to be expected for higher values of h/5j. 

The static pressures and the drag coefficients corresponding to these static pressures 
which result from turning the free-stream flow through the separation angle 6 ^ are also 
shown m figure 2. This makes it possible to determine the relative contribution of the 
static and dynamic pressures to the drag coefficient and dividing-streamline stagnation pressure. 
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The results indicate that the contribution of the dynamic pressure is low for low Mach 
numbers and low values of h/5j, and high for high Mach numbers and high values of h/5^. 

Effects of changes m 6^ and a.- In order to apply the theory, two empmcal 
factors, the separation angle and the spreading factor a, are used. These are based on 

correlations of questionable accuracy; consequently, figures 3 and 4 have been prepared to 
indicate the sensitivity of the present method to changes in these two parameters. Calcula- 
tions are presented in figure 3 for 9^ -t- 1° and 6^ - 1° and in figure 4 for 0 9a and 
l.lo, where the values of 9^ and a are adopted from reference 1. 

A decrease of 2° in 9^ (fig. 3) results m a significant decrease in the dividing- 
streamhne stagnation pressure and in the drag coefficient. The effect is largest for the 
highest step-height ratios and free-stream Mach numbers 

An increase of approximately 20 percent in the value of o (fig. 4) decreases appre- 
ciably the values of the dividing-streamlme stagnation pressure and drag coefficient. The 
strongest effect occurs at the lowest step-height ratios and highest free-stream Mach numbers. 

A study of more recent correlations for o (ref. 3) indicates that a probably should be 
a function of h^Sj, Rgj) ^nd Mach number and this could have some influence on the 
shape of the drag curves for h^5j > 5. 

Effects of changes in velocity profile .- The effects of the changes in the power of the 
velocity profile from 1/7 to 1/9 to 1/11 are illustrated in figure 5. As expected there is 
only, a small effect of the change of the velocity profile at the higher values of h/Spbut 
the effect becomes stronger as h/5j decreases For Pod/Pp greatest change, an 
increase with increase m power index, occurs at the highest Mach number, for Cj^, the 
largest change, a decrease with increase m power index, occurs at the lowest Mach number 
The maximum vanation in p^^^Pj and Cj-j is less than 8 percent 

Pressure distributions .- The pressure integral used to calculate the step-drag coefficient 
implies a certain pressure distribution on the step face. This pressure distribution is made 
up of two components, the static pressure due to turning the stream flow through the 
separation angle 9^ and the pressure derived by stagnating the flow in the shear layer as 
well as reverse-flow region. Some typical implied pressure distributions are shown in figure 6 
for a range of ratios of step height to boundary-layer thickness for Mach numbers of 3 and 6. 
For reasons that will be discussed later, the assumption was made that there were no normal 
shocks in the flow even if the flow was somewhere supersonic. For a particular step-height 
ratio the pressure distribution region is divided into two parts. One is in the shear-layer 
region and another is in the reverse-flow region. In the shear-layer region, the step-face 
pressure decreases from the dividing-streamline stagnation or peak pressure nearly exponentially 
toward the static pressure of separation. In the reverse-flow region, the step-face pressure is 
uniform everywhere because of the assumption of uniform velocity and uniform static 
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pressure behind the separation shock. There is a discontinuity in the flow at the junction 
of the shear layer and reverse-flow region. The results also indicate that the step-face pres- 
sures are of higher magnitude for higher Mach numbers with the dynamic pressure component 
increasing in importance. 

As h/5j decreases, the dividing-streamline stagnation pressure decreases, but the shear- 
layer height on the step face increases. Thus the reverse-flow height on the step face also 
decreases and results in higher dynamic pressures in this region. For very small values 
of h/6j, the pressures m the reverse-flow portion exceed the peak pressure at the outer 
corner. (See fig. 6.) This is physically impossible and explains the rapid rise of the curve 
of drag coefficient with decreasing h/Sj- ^*8- 2(b).) Use of energy conservation 

would probably yield an improved estimate of the pressure in the reverse-flow region but 
with an increase in complexity. 


Variable 

Pressure distribution .- The previous problem is bypassed by fixing the junction of the 
shear layer and the reverse-flow region at the center of the step face. For this condition, 
the IS not kept constant, but is allowed to vary so as to satisfy the previous condi- 

tion. The resulting pressure distributions calculated by this method are shown in figure 7 
for Mach numbers of 3 and 6 for the same range of h/S^ as in figure 6. In the shear- 
layer region, the pressure distributions are almost the same as those in figure 6, but in the 
reverse-flow region, the pressure distnbutions are essentially independent of h/5|- This 
characteristic tends strongly to delay the breakdown of the numerical calculations for low 
step-height ratios, but the theory must still ultimately deteriorate as h/8^ decreases 
because the pressure integral (eq. (3)) is singular for h^5j = 0. 

Reattachment pressures and drag coefficients .- Figures 8(a) and 8(b) show comparisons 
of the dividing-streamline pressures and drag coefficients calculated by both the variable and 
constant methods. For both methods, the pressure ratio Pod/^i computed 

by means of the isentropic equation (6). For the case of Pod/Pp signif- 

icant difference between the results of the two methods for h/5j > 1,0. However, 
for h^8j < 1.0, the vanable method generates Pq d/Pi values that are higher than 

those calculated by the constant method. The variable method results in some- 

what higher Cj) values at the higher values of h/5j and lower values at lower h/5j 
values, cumulating in a curve that shows a decreasing drag with decreasing h/5j and is 
more in agreement with experimental trends Also, the tendency of the theory to break 
down ^as evidenced by the upturning of the drag curves at the lower values of is 

delayed to at least an order of magnitude decrease in step-height ratio and suggests the 
possibility of extending the usefulness of the theory to these lower values of h/5j. 
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Reattachment pressure deriv ed from pressu re j^ntegral.- The use of Crocco’s integral or 
relationship between velocities and densities with the static pressure and stagnation tempera- 
ture held constant yields the relationsliip 


^a 


1 - c: 


2^2 


c:0 


( 10 ) 


The use of the isentropic relationship in equation ( 6 ) (for shockless flow) 


Pi 



+ 




Pi 


cannot satisfy both these requirements simultaneously. For a constant-stagnation-temperature 
adiabatic process, the density relationship can be shown to be (ref. 4) 


^a 



( 11 ) 


In equations (10) and (11), 02 /P 3 increases as 0 increases as required in the shear 

layer, but the rates of increase are completely different due to the different power index. 
Thus, Pod/Pi calculated by equation ( 6 ) is not compatible with the problem. For 
this reason the reattachment pressures were also calculated directly from equation (7) 

^with 0 = 0j^ which was denved from Uebelhack’s pressure integral. The results are 
presented m figure 9 and are compared with present calculations made by Uebelhack’s 
method (eq. ( 6 )). 

The dividing-streamline stagnation pressure calculated by Uebelhack’s method is not in 
agreement with that calculated by the integral equation. There is considerable difference 
in values at low Mach numbers and low values of h^Sj. The integral method gives pres- 
sure ratios p^ ^^p ^ that are higher than those calculated by the Uebelhack’s method for 
the lower Mach numbers and the ratios increase more slowly with the increase m Mach 
number. With the exception of M, = 6 , the integral method gives higher values of pres- 
sure than Uebelhack’s method for the range of h/5j considered. 

Mach number on dividing streamline .- The Mach number on the dividing streamline 
is of interest since it indicates the type of equations required to determine the conditions 
within the separated-flow region for more exact methods Consequently, figure 10 shows 
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the dividing-streamline Mach number Mj at the outer corner of the step, both as a func- 
tion of Mj ^ and of h/5j. The figures indicate that increases as ^ and h/S^ 

increase. For Mj ^ > 3, is supersonic for h/6j > 4. For Mj ^ = 5, is 

supersonic for h/6j > 1. These plots indicate that for many supersonic free-stream Mach 
numbers and step-height ratios, one should consider the possibility of a normal shock near 
the step face if accurate pressure distributions are required. 

Effect of normal shock. - The effect of a normal shock on the flow is illustrated in 
figure 1 1 . The most severe practical case for the present range of calculations is considered 
^Mj g 6 and h/6j = 10^. 

The effect of the normal shock is significant on the pressures near the dividing 
streamline when the flow is supersonic, but the effect on the integrated drag coefficient 
IS only about 3 percent if the same total pressure loss is assumed in the reverse-flow region 
as m the shear layer. This error in the drag coefficient is so small that it was not con- 
sidered necessary to include the normal shock effect in the drag calculations in this report. 
Due to the normal shock, the static pressure at the step face is drastically changed and is 
no longer constant. Consequently, the assumption of constant static pressure everywhere 
behind the oblique shock is not valid with the normal shock present and any attempt to 
include this variable static pressure m the theory greatly increases the complexity of the 
calculations. 

The largest error ACq/Cq resulting from negligence of the normal shock effect 
on Cp occurs at the higher Mach number investigated (Mj ^ = 6^ and is shown as a 
function of h/6j in figure 12. The ratio ACp/Cp is the difference in the integrated 
drag coefficient with and without a normal shock existing whenever the flow is supersonic. 
The curve indicates that ACp^Cp increases rapidly from 0.001 at h/5j = 1 to an 
asymptotic value of about 0.04 for h/5j > 100 Thus, the effect of the normal shock 
on Cp can be neglected, for the range of conditions covered in this report. 

Comparison With Experiment 

Pressure distributions .- Because of the great number of simplifications and assumptions 
involved in the present analysis, the accuracy and usefulness of the theory can only be 
established by comparison with experiment Confidence in the theory can be enhanced 
only if it can be shown that the pressure distributions obtained by the theory are in 
reasonable agreement with expenmental ones. Comparisons are presented in figure 13. Inas- 
much as the experimental data are known not to have a constant pressure distribution on 
the step face in the reverse-flow region, the theory is modified at this point with the 
assumption that the pressures in the reverse-flow region are a mirror image of those in the 
shear-layer portion of the step face with the focal or turning point located at half the step 
height. For this special case, the requirements of conservation of mass and energy are 
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automatically satisfied and the momentums in the two regions are equal although in opposite 
directions. Consequently, this special case meets much more stringent requirements than the 
general theory where only conservation of mass is involved. 

In general, the theoretical pressure distnbutions have a reasonable similarity to the 
experimental ones and are of the correct order of magnitude. The comparison suggests, 
however, that the assumed separation angles 0g are too small at all Mach numbers. An 
increase in 6 ^ will also tend slightly to improve the agreement m the shapes of the 
theoretical and expenmental curves in the shear-flow region. The expenmental data indicate 
that the pressures in the reverse-flow region are not as high as those in the shear-flow area. 
This trend suggests that the effective static pressure, due to 0^, is not constant but decreases 
with distance from the upper corner of the step. The minimum value in the expenmental 
pressures appears to occur in the general area of y/h == 0.4 rather than the value of 0.5 
assumed in the theory. Note, finally, that many of the experimental pressure distributions 
are not truly adequate to determine reliable values of either the dividing-streamline stagnation 
pressure or step-face drag coefficient. 

Reattachment pressures and drag .- A comparison of the theoretical and experimental 
values of Pod/P i corresponding to the pressure distributions of figure 13 and 

additional available experimental data are presented in figure 14. Experimental data are 
taken from references 1 and 5 to 10. Table I presents the relationship between the symbols 
and basic parameters for the experimental data along with the sources from which the data 
were obtained. In some cases the extrapolations of some of the experimental data were so 
large that the p^ ^ and Cjj values that were determined were considered to be too 
unreliable to justify plotting. In some instances the experimental Cp values determined 
from experimental pressures do not match those determined by Uebelhack. In particular, 
the data of Sterrett and Barber (ref. 5) as presented by Uebelhack appear to be m error. 
Furthermore, some of Uebelhack’s data could not be plotted as a function of h/6j because 
insufficient information was given for identification. 

The comparison indicates that m general the theoretical reattachment pressures are in 
good agreement with expenment as far as magnitudes of values and trends with Mach 
number and step-height ratio are concerned. There is some scatter m the expenmental data 
with the data of Uebelhack appearing to be consistently lower than the data from the other 
sources. The data of Sterrett and Barber also appear to indicate a much more rapid rise 
in Po,d/Pi than predicted by the theory for the higher experimental step- 

height ratio. More data are required to establish whether this is a true trend or merely 
an expenmental abenation. 

The drag results, on the other hand, indicate that the theory tends to generally under- 
predict the experimental values but predicts the correct variation with Mach number. The 
theory also probably predicts the correct variation with h/5j; however, the expenmental 
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data are too meager and the scatter too large to verify this fact satisfactorily. The inability 
of the theory to predict the correct magnitude of the drag coefficient is ascribed to the use 
of too low values of the assumed separation angle over most of the Mach number range. 

In view of the erroneous density relations used to determine Pod/Pp ^^e improper 
interpretation of Sterrett and Barber’s results, and the incorrect calculated trend with h/5j, 
the agreement between theory and experiment presented by Uebelhack appears to be some- 
what fortuitous. 


Empincal Extension of Theory to Small Step-Height Ratios 

The theory was developed for cases where h/6j > 2, however, there is a great need 
to develop drag-prediction methods for h/5j < 2. Consequently, an investigation was made 
to determine the requirements for developing such an extension using some unpublished 
experimental data together with data from reference 6 as a guide. 

The unpublished data (see ref. 1 1 for some limited results) indicated that the step-face 
pressures and drags decreased to 0 as h/6j approached 0. Furthermore, the data indicated 
that both the static pressures due to the separation angle 0g and the dynamic pressure m 
the shear layer had to decrease with decreasing h/5j. In order to achieve this trend, a 
decrease in 6^ had to be coupled with an increase in the shear-spreading parameter a. 
Obviously, a large number of combinations of 0g and a could provide the desired 
results. In order to keep the calculations orderly, the assumption was made that there was 
a functional relationship between 0g and a and that this relationship could be expressed 
as 


od 


s 



( 12 ) 


The relationship between 6g and h/6j was assumed to be given by 


®s = 




(13) 


This expression provided the best fit to the experimental data of the 8 types of expressions 
investigated. Equation (13) allows 0g to decrease very slowly at first from the reference 
value at some chosen value of h/6j, then decrease much more rapidly at the lower step- 
height ratios, and end with a slow decrease in 0g toward 0 as h/Sj tends toward 0. 
Calculations were then made at Mj ^ = 1.61 to determine the boundary conditions that 
would fit the empirical curve to the experimental data. These boundary conditions are 
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Upper boundary conditions. 


— = 7 0 = 12.25° 

Lower boundary conditions' 


30s 

Wi) 


(14) 


90 

— = 0.004 6, = 0.004" . . = 0.004 or 0.0 (15) 

S, ' 4(l>/«l) 

The same boundary conditions except for upper 0„ (for which 13.25° was used) were 
then used to calculate the theoretical curve at M, = 2.2. The results of all these calcu- 

lations are shown in figure 15 where both Pod/P 1 plotted as functions 

of hySj. The experimental data are shown in the form of shaded bands. The widths of 
these bands are pnmanly due to Reynolds number effects and not scatter. These Reynolds 
number effects can be readily incorporated into the theory; however, no attempt was made 
to do so in the present investigation. These experimental Reynolds number effects are 
considerably larger than the theoretical effects calculated for changes in the index of the 
power-law velocity profile, and hence suggest then an additional parameter not considered 
herein may be involved. 

Figure 15 shows that, in general, theory is in good agreement with experiment at 
both Mach numbers. Note that this agreement extends down to h^5^ in the order 

of 0.001 for the drag coefficients despite the fact that the pressure integral (eq. (3)) is 

singular for h/5 j = 0. The experimental data for Pod/^i were not plotted to such 
low values of h^5 j because the models had become so small m height that, in the 
extreme case, only one or two orifices were available to represent the pressure distributions, 
and pressure peaks could not be isolated. Note also that the experimental data never fall 
below the static pressure curves and thus appear to confirm a reasonable division between 
static (0j) and dynamic pressure (a) effects in the extended theory 

The value of h^6^ = 7 required for the upper boundary condition m matching the 
theory with experiment may appear to be too high and the possibility exists that the 
function chosen to relate 0^ to hy'Sj may be too elementary and the approach of the 
extended theory to the basic theoretical curves may occur more sharply. Attempts to resolve 
the problem, however, uncovered just as many reasons why the value may be valid as why 
the value is too high. 

Although the theory appears to be in good agreement with experimental data for h/6j 
ratios as low as 0.001, recommendations are made that this theory not be applied to h/5j 


17 



values less than about 0.05 because of the increasing Reynolds number effects, which are not 
incorporated into the empirical extension of Uebelhack’s basic theory. 

Finally, a comparison is made of some pressure distributions predicted by the extended 
theory with some unpublished experimental results obtained at Mach numbers of 1.61 
and 2.20 (fig. 16). An effort was made to pick expenmental data that was representative 
of the central regions of the shaded areas in figure 15 and also reasonably close to the rela- 
tive location of the theoretical lines in the shaded areas. For the results at a Mach number 
of 2.20, the Rgj = 0.742 X 10^ data correspond more closely to the center regions of 
the shaded areas and the R 51 - 0.524 X 10^ data relate more closely to the areas next 
to the theoretical curve in the h^5j range of interest. As before, the theoretical reverse- 
flow pressures are assumed to be mirror images of those m the shear-layer region. Note that 
the theoretical and experimental step-height ratios do not coincide exactly. 

In general, the agreement between theory and experiment is good. A significant feature 
indicated by the experimental results is that the reattachment point (or the dividing stream- 
line) does not occur at the upper corner of the step but at some distance below it and that 
this distance increases as h/5j decreases. Behrens (ref. 7) also notes this effect. What 
this means is that the shear-layer flow above the reattachment point passes over the top of 
the step and does not enter the separated-flow region. Consequently, for a more correct 
companson with experiment the theoretical region should be foreshortened to include only 
the distance from the wall to this peak in pressure. Such a foreshortening will substantially 
improve the agreement between the theoretical and experimental pressure distributions and 
move the theoretical pomt of minimum pressure closer to the y/h region of 0.40 to 0.45 
where the expenmental minimums generally occur. Of importance is the fact that most 
experimental results do not contain sufficient detail of the pressure distributions near the 
upper comer, making it difficult to accurately estimate reattachment pressures or face-drag 
coefficients. Finally, the agreement between theory and experiment would not have been 
quite as good if no effort had been made to match the data approximately for Reynolds 
number effects. 

Another point of interest revealed by this analysis is that as h/5^ decreases the com- 
ponent of pressure contributed by the dynamic pressure decreases much more rapidly than 
the component contributed by the static pressure rise across the separation shock. Finally, 
mention should be made of the fact that the static-pressure measurements of Behrens (ref. 7) 
on the wall ahead of the step and also outside the shear layer indicate that the linear shear 
layer probably does not exist for h^S^ < 3. Thus, it is quite surprising that the extended 
theory appears to perform so well. 
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SUMMARY OF RESULTS 


A theoretical investigation has been made of the pressure distributions on and drag 
charactenstics of the faces of forward-facing steps immersed m turbulent boundary layers 
at supersonic speeds Uebelhack’s theory was used as a basis for the calculations but the 
theory was modified and extended, and a different method was used to obtain the solutions. 
The results of this investigation are summarized as follows 

(1) Decreases in the assumed separation angle and increases in the assumed shear-layer 
spreading parameter result in decreased pressure levels and drag forces, with the separation 
angle affecting the effective static pressure and the spreading parameter the dynamic pressure 
m the separated-flow region. 

(2) Increases in the index of the assumed power-law velocity profiles caused small 
increases in the dividing streamline or reattachment pressures at the smaller ratios of step 
height to boundary-layer thickness but had no significant effect on the drag. 

(3) Modification of the theory so that the shear layer always terminates at the center 
height of the step face and replacement of the rectangular pressure distnbution of the reverse 
flow with the mirror image of the shear flow resulted m increased accuracy m the calcula- 
tions and in reasonably good agreement between theoretical and experimental face pressure 
distributions for moderate ratios of step height to boundary-layer thickness (from 2 to 7). 

(4) By coupling a decrease in assumed separation angle with an increase in the shear- 
layer spreading parameter, the modified theory can be extended to give good agreement 
between theory and experiment to ratios of step height to boundary-layer thickness decreas- 
ing to values as low as 0.05 at free-stream Mach numbers of 1.61 and 2.20. 

(5) There is some disagreement between the results obtained in this investigation and 
those of Uebelhack which apparently cannot be ascribed to different methods of computation. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
August 26, 1975 


19 



APPENDIX A 


THEORY 

The primary assumptions m Uebelhack’s shear-flow model for boundary-layer separation 
ahead of a forward-facing step have been mentioned in the section entitled “Theoretical 
Approach.” For details on the secondary assumptions and general development, see refer- 
ence 1. ■ The main objective herein is to point out some deviations from Uebelhack’s 
approach that were utilized in the present paper and to correct the typographical errors 
in some equations in reference 1. 


Separation Region 

The boundary-layer momentum thickness ^2 behind the oblique shock at the separa- 
tion point is an important parameter combining mass and momentum flux for the free-mixing 
region and is calculated as follows (see fig. 17): 

The conservation of mass in the boundary layer across the oblique shock is given by 

lp,Ui dy dy (Al) 


The momentum flux across the oblique shock in the x-direction is given by 
Pl^l ^^l^l ^^2^1 ’’’ 


The oblique-shock theory is then applied to calculate 


_ sin 


^2,e^2,e 


(<5 - «s) 


sm /3 


and 


<‘l,e“l,e , ^1^ , s.n(g-9,)cos(g-«,) 

2 , ,2 sin B cos (3 

^2,e^2,e P2^2,e 


(A2) 


(A3) 


(A4) 
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Introduce the boundary-layer displacement and momentum thicknesses defined by; 


6 . . p(i . 

•'O ' 


dy 


(A5) 


e = /■* -a_A . 

•'0 ^e^e ' 


(A6) 


Use of equations (A3) to (A6) in equations (Al) and (A2) allows the calculation of the 
boundary-layer momentum thickness as follows’ 


= 

2 cos 6c 


(5, - - «,) 


sin (|3 - cos (|3 - 


sin j3 


cos j3 


(«i - «i) 


sin (/3 - 0„) 


’1 - I ^ cos 0„ + 5, 

1 V sin j3 s 1 


Pi - P2 




(A7) 


Free-Shear Layer 

The shear-layer velocity profile in the mixing region is approximated by the error 
function profile of the form 


U2 U2 j 

— = = <f> = + erf (tj)] 

“a “2,e 2 


(A8) 


where 


erf (t?) = 


n/tt Jq 


2 

e - d? 


(A9) 


_ oy 

V = — 

X 


(AlO) 


where ^ is a dummy variable. 

The spreading factor a is not known exactly for compressible fluids; Korst and 
Chow (ref., 2) suggest the following correlation 


(7 = 12 + 2 76M, 


2,e 


(All) 
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For the constant stagnation temperature and constant static-pressure flow with unit 
Prandtl number, Crocco’s integral relationship is applicable. In this case, the density ratio is 
given by 



As the shear layer grows in the downstream direction from the virtual origin, the 
origin of the coordinate system, which is anchored to the center of the velocity profile, is 
displaced toward the lower velocity. The coordinate shift yj^ is determined from the 
continuity and momentum equations using the control volume indicated in figure 18. The 
continuity equation is 


-Pb^byrb + Pa^ayra P\,\ 



pu dy 


(A 13) 


The momentum equation is 


-^b^^rb 


2 

p u y 
a*^ra 


r 

•'n 


'=b''b“b 



dy 


(A 14) 


In these equations yj.^^ and yj.^ are arbitrary large negative and positive values, 
respectively, such that 


u 



u 


2,e 


and 


^(yrb) ^ ^b 

and the integral value does not change significantly for further increases in the limits. After 
combination of equations (A13) and (A14), normalization with P^u^, and with the setting 
of U|j = 0, the result becomes 

ym = yra - f- dy (^15) 

•'yrb 
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By use of the definition of the equations (A 10) and (A 12), equation (A 15) becomes m 
nondimensional form 


’^m ^ra 



(' - c\y 
1 - 


drj 


(A16) 


The jet boundary streamline rjy which separates the original jet from the entrained 
part of the flow, for this case of rectangular initial velocity distribution can be calculated 
from the continuity equation (see fig. 18) 




I pu dy 

‘'vi 


(A17) 


After partial integration and conversion of the above equation to dimensionless form, it 
becomes 


^ra ^m 




dr/ 


(A18) 


With the use of equation (A16), equation (A18) is further rearranged to calculate the jet 
boundary streamline 



(A 19) 


In the absence of a boundary layer, the tjj streamline is identical to the dividing 
streamline rj^. When an initial boundary layer exists, the dividing streamline must be 
displaced to lower velocities in the shear layer to allow for the momentum losses in the 
boundary layer. This effect can be created by injecting a mass flow of zero x-momentum 
into the shear flow. (See fig. 19.) From the continuity equation, the mass m^jj between 
the jet boundary streamline tJj and the dividing streamline rjj is given in dimensionless 
form by 
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^ a, - d. 

w 1 - cy > - cy 


(A20) 


Note that the location of the dividing streamline relative to the boundary streamline is 
defined only by the continuity equation; the momentum equation is not taken into account. 
This approximation is reasonable for the small mass bleeds involved. The dividing streamline 
is obtained from the conservation of mass and momentum in the control volume of 
figure 19. It IS found to be 



^ nr. (■ - cl)o (■ - *) 

\h 1 - 


dr? 


(A21) 


where 02 the boundary-layer momentum thickness. 

Due to the injection of the air in the shear layer, the shear layer is thickened. Hence 
the virtual origin is moved ahead of the actual separation point along the center line of the 
shear layer (fig. 20). Thus the length of the shear layer is increased in the system by the 
amount (from that without boundary layer), and the new length of the shear layer 

X is given by 

X = X* -I- Xq (A22) 


The analysis for calculation of the length x^ is shown below. 

The term a02/x* determines tj^ in momentum equation (A21) and the term 
am^jj^p^u^x* determines in continuity equation (A20), so that simultaneous solution 

of the two equations gives 

”"dj = Pa^a^2 ^^^3) 

The mass m^j can be considered to be injected and entrained between x = 0 and 
X = Xq. Hence 

"^dj = Pb^b dx = PbVbXo (A24) 

*'o 
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From equations (A23) and (A24), Xq can be related to $2 


^o _ ^a^a _ 1 ^a 

^2 1 - '"b 

Also, for the control volume of figure 19, the continuity equation is used to determine 
the injected mass below the jet boundary streamline i?j Thus, another expression for Xq 
is obtamed 


X 


o 


^2 


a 


•>rb ‘ - 


drj 


(A26) 


For the new arrangement, a new dimensionless lateral distance r? is defined as yo/x 
and the same theory as for the shear layer without injection is utilized. The dividing- 
streamline location in the barred r? coordinate system is obtained from equation (A21) as 


pi (l - ^ + pa 

’'^rb ^ ^ T^h 


1 - 0^02 


(A27) 


where the integral limits 7 ?j,g and are arbitrary numbers sufficiently large so as not 

to affect the value of the integral significantly. The integral value, hence, will not change 
by changing the limits rjj.^ and to Tjj.g and For actual calculations, the 

integral on the right side of equation (A27) is replaced by the term a02^XQ which can 
be calculated in the original unbarred coordinate system from the following equation derived 
from equations (A26) and (A 19) 


riin, (■ - Cj)»(l - ») ^ p (■ - ^ ni (■ - Cl)0 

1 - c|02 •'brb 1 - 


002 


(A28) 
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Reattachment Region 

The point of reattachment of the dividing streamline is assumed fixed at the outer 

comer of the step. Because of the adjustment of the shear layer for the initial boundary- 

layer thickness, the dividing streamline no longer hits the corner of the step as originally 
assumed. Consequently, the new location of has to be found by iteration until the 

dividing streamline meets the corner of the step within the desired accuracy. 

The impingement point of the dividing streamline on the step face is usually initially 

located below the corner by the distance At? (see fig. 21). This distance can be reduced 
to zero by the equation (derived from the geometry of fig. 20) 


Ax = At? — tan (A29) 

o ^ 

where 

Ap = Pm - ^d 


for the first iteration and 

Ap - P(j(old) ~ ^d(new) 

for subsequent iterations. Similarly, the longitudinal distance x is computed for each 
iteration from 

^new = ^old + 


A new value of pj that is closer to but still not at the step corner is, of course, calculated 
from equation (A28) each time a new value of x becomes available. The iteration process 
is continued until Ap becomes insignificantly small m relation to the step height. 

The Mach number on the dividing streamline is given by 





1 - 




(A31) 


Uebelhack calculated the reattachment stagnation pressure p^ for this streamline by the 
isentropic one-dimensional flow equation, where < 1 
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~ 1 ^ — (A32) 

2 d/ 

and by Rayleigh’s formula when M^j > 1. The method of calculation of Po,d/Pi 
subsonic is not correct, because the flow is stagnated by the isentropic process at the 

step face whereas the assumption has already been made that the flow in the shear layer 
has constant static pressure and constant total temperature. For Mj > 1 there is insuf- 
ficient information given to make a comment at this point. For the correct method of 
calculating the desired pressure ratio some additional denvations are required, hence the 
method will be discussed later. 


Pq,.! _ Po,d P2 - (i I 
P, P, P. ' 


Reverse-Flow Region 

The shear-layer flow and the reverse-flow regions are shown in figure 22. The height 
of the shear layer below the dividing streamline is AB and AE is the vertical total 
distance for the reverse-flow region. The distance AB is given by 


AB = 

and the distance AE is given by 


(A33) 


AE = h - AB cos = h - - '^rb)'^ ^s (A34) 

The mass flows through AB and AE are assumed equal. With this assumption, the 
reverse-flow velocity can be found from (with the reverse-flow velocity assumed uniform) 


^a 



£ 


T?d (l - c2)0 


^rb 


- cy 



1 

(^d - ^rb) ^s 


(A35) 


With the velocities and densities now known in the region near the step face, the , 
momentum equation can be applied to the control volume shown m figure 19. Lines AB 
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and AE are assumed close to the step face so that the shear term can be neglected. 
After some manipulation and simplification, the stagnation pressure integral becomes 



P - P2 


dy 



n 2 
2 

Pa^a 


dy cos + 



2 

Pb% 


2 

Pa“a 


dy 


(A36) 


The expressions for the density ratios ^ 2/^3 ^b/P^ substituted from equa- 
tion (A12) and substituted from equation (A35) into the equation (A36) to 

yield 




2 


dr} cos 


+ 






X 


' 

1 

(^d - ^rb) cos 

I 


(A37) 


The pressure integral is related to the free-stream static pressure pj existing ahead of the 
separation shock by 



iy = 

Pj h 



P-^2 dy ^ 

P 2 ^ 



The drag coefficient on the step face is calculated from 


(A38) 


C 


D 



_P. 

Pi 

-^tM 

2 



2 

l,e 


(A39) 
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The correct dividing-streamline pressure is calculated by the use of nght first term 
of equation (A37) after the removal of the integral sign which yields 



(A40) 


Pp,d 

Pi 



(A41) 


Moreover, the stagnation pressure distribution anywhere through the shear layer below divid- 
ing streamline (i.e., along AB) can be calculated from equations (A40) and (A41) by 
changing to (p and p^ ^ to p. The Mach number along the shear-layer line AB 

IS calculated by 




2 2 
Zi 


1 - C 


2.2 


(A42) 


When the free-stream Mach number and h/5j are sufficiently large, a part of the shear 
layer in AB is supersonic. For this case the flow was decelerated through a normal 
shock with the stagnation pressure behind the shock being found by using the well-known 
Rankine-Hugoniot relationships (ref. 4). 
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CALCULATION PROCEDURES 

The calculation begins with the specification of 5j, h, and n. With Mj 

as the entering variable, the separation angle is found from figure 1 in reference 1 . 

From the shock angle /?, the Mach number behind the shock, and the static pressure 
rise P 2 /P 1 found from tables such as those in reference 4. Two other required 

quantities and can be found in terms of the ratios ^l/^l 

tables in reference 12 as functions of M, and the power of the velocity profile. 

1 

Finally, the spreading factor o is calculated from equation (All) in appendix A and the 
Crocco number C^ is determined with the aid of equation (51) of reference 4 utilizing 
the Mach number behind the separation shock as the reference condition. 

With all the required basic quantities now known, the next step is to calculate 02 
by equation (A7), by equation (A16), and by equation (A 19) In equa- 

tion (A 19), rjj IS a variable upper limit m the integral in the left-hand side of the equa- 
tion whose value must be found to satisfy the values of the two integrals on the right- 
hand side. For the present calculations, the limits tJj.^ and were chosen as 3.5 

and -1.5, respectively, and calculations indicated that the values of integrals involving 
these limits were not affected significantly by changes in these limits provided that the 
changes were m the range or > 1.5. With the use of rjj, the value 

of Xq is calculated from equation (A26) and thence the value of x from equation (A22) 
This value of x is then introduced into equation (A27) to determine the initial value 
of Tjj. The final value of 17 ^ is found by iteration with the use of equations (A29) 
and (A30) to find a new value of x which is reintroduced into equation (A27) to find 
a new value of The iteration is continued until is established within the limit 

of desired accuracy. Finally, with r}^ known, the pressure integral is computed with 
the aid of equations (A37) and (A38) and the drag coefficient from equation (A39). Note 
that this procedure does not allow for the existence of any normal shock within the 
separated-flow region. 

The procedures for calculating additional items of interest such as M^j, Po,d/Pl> 
and pressure distribution are quite straightforward and are not detailed here. 

The calculation time for a typical run which might include 8 step-height ratios at 
5 Mach numbers for a total of 40 points at moderate step-ratios (0.5 to 10) was on the 
order of 200 sec CPU on the Control Data 6400 computer system. 
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A Separation region 
B Free-mixing region 
C Reattachment region 
D Reverse-flow region 


Figure 1.- Basic flow model and subregions. 
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Figure 8.- Comparison of results computed by constant and 
variable ^ methods. 
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Figure 13.- Comparison of theoretical and experimental pressure distributions. Variable 
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